Origin of the spatial variation of the pairing gap in Bi-based high temperature 

cuprate superconductors 
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Recently, scanning tunneling microscopy on Bi-2212 cuprate superconductor has revealed a spatial 
variation of the energy gap that is directly correlated with a modulation of the apical oxygen position. 
We identify two mechanisms by which out-of-plane oxygens can modulate the pairing interaction 
within the Cu02 layer: a covalency between the x^-y^ band and apical p-orbital, and a screening 
of correlation U by apical oxygen polarization. Both effects strongly depend on the apical oxygen 
position and their cooperative action explains the experiment. 

PACS numbers: 74.72.Hs, 74.62. Bf, 74.20.-z, 75.30.Et 
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Various types of high-Tc cuprates have been discov- 
ered for the last few decades. As a function of density of 
charge carriers doped into the Cu02 planes, a supercon- 
ducting critical temperature Tc shows in general a max- 
imum which varies strongly from one family of cuprates 
to another. For hole-doped cuprates, a correlation be- 
tween the maximum Tc and the energy-level separation 
of in-plane oxygens (Op) and apical ones (Oa) has been 
noticed early on [J, [3] . Hence, a role of apical oxygen on 
Tc is of considerable interest [l|, S i, i, 0, @, • 

In Bi2Sr2CaCu208 (Bi-2212), a mismatch between the 
rock-salt Bi02 layers and the Cu02 planes causes an ex- 
tra modulation of the crystal structure with a period 
about 26 A. In such a " supermodulated" lattice, the dis- 
tance d from the Cu02 layer to the apical oxygen Oa is 
periodically varied within the range of ^ ±6% [8[ . Given 
that bonds within the Cu02 plane itself are much less 
affected by supermodulation, Bi-2212 material provides 
a unique opportunity to study the impact of apical oxy- 
gen on superconductivity, by monitoring local electronic 
properties as a function of d that varies spatially within 
a supermodulation period. This is precisely what is done 
in the recent scanning tunneling microscopy (STM) ex- 
periment by Slezak et al. [l3| (see also Ref. [ll|) . 

Slezak et al. measured the gap A in a local density of 
electronic states, and found nearly 10% spatial variation 
of A with the same periodicity as supermodulation. They 
have emphasized that the gap variation is anticorrelated 
with the Cu-Oa distance variation, i5A oc —5d, i.e., the 
gap increases when Oa gets closer to the Cu02 plane and 
vice versa. 

This remarkable observation has already been ad- 
dressed in several papers by introducing variations of 
coupling constants on a phenomenolog ical level [13, El- 
However, the underlying microscopic mechanism that 
links the strength of pairing interactions within Cu02 
planes with the position of out-of-plane Oa remains elu- 
sive. In this Letter, we discuss the physical origin of high 



sensitivity of the pairing gap to the Cu-Oa distance and 
explain the anticorrelation effect (5A oc —5d observed. 

Quite in general, the structural shift of apical oxygen 
may influence the energy gap either via the hopping pa- 
rameters (hence density of states on the Fermi level), or 
via the strength of the pairing interaction. In particular, 
there is a well-known relation between the next-nearest- 
neighbor hopping t' and Tc, based on the band structure 
calculations by Pavarini et al. One has to notice, 

however, that this observation concerns a comparison be- 
tween different families of cuprates with different lattice 
structure. In fact, Pavarini et al. predicted that the vari- 
ation of the Cu-Oa distance within a qiven compound 



hardly affects the hopping parameters 1J|. Therefore, 
we focus here on possible electronic mechanisms by which 
apical oxygens may affect the strength of the pairing po- 
tential within the x'^-y'^ band [ist . 

We find two different ways how the apical oxygen may 
enter the game. First, a hybridization of "useful" x'^-y^ 
band with the " pairing- inert" orbitals of apical oxygens 
reduces the pairing interaction. Such a destructive effect 
of covalency is controlled by a relative energy separation 
between the orbital levels that depends on Cu-Oa dis- 
tance. We illustrate this by an explicit calculation of the 
Madelung potential as functions of d. Second, we show 
that the superexchange interaction J, which is believed 
to be essential for magnetic correlations and possibly for 
superconductivity, is very sensitive to the apical Oa po- 
sition. This is due to the high polarizability of 0^~ an- 
ion which has the effect of screening and reducing the 
energy U needed to move an electron from one ion to an- 
other [il,[l3. InB i-2212, the screening effect and hence 



the strength of magnetic correlations J (x 1/U are spa- 
tially modulated because the closer the apical Oa is, the 
stronger the screening is. Remarkably, we find that the 
above two effects, covalency and screening, both favor an 
anti-phase relation between A and d variations. 

Covalency.- We address this effect in terms of the fol- 



lowing Hamiltonian: 

ka k 

+ Xl'^AaLafe'T + ^J^fc(aLcfe^ + (1) 

ka k{T 

Here, the first term corresponds to tlie in-plane pda-band 
of x^-y^ symmetry made of Cu and Op orbitals, and a 
conventional form of dispersion 4] = — 2t(cos + 
cos ky) + At' cos kx cos ky — 2t"{cos2kx + cos2ky) — fi, 
where fj, is the chemical potential, is adopted. Second 
term shows that this band is supposed to host a su- 
perconductivity of c?-wave symmetry, with the gap func- 
tion 7^ = (cosfca; — cos ky)/ 2 and the gap magnitude 
A = 5 X]/c 7fc ('^-fciCfe-f }, determined by the strength g of 
a pairing potential (whose origin is not specified). 

Third term in Eq. ^ represents holes on the 2pz or- 
bital of apical Oa with energy ea. While a complete 
model may consider the rich internal structure of the "ax- 
ial" orbital by Pavarini et al. including Sd^z^-r^ and 4s 
states of Cu [4|, we consider here a minimal model that 
captures the essential effects of the axial orbitals. Fi- 
nally, the last term accounts for a covalent mixing of the 
x^-y^ and pz orbitals. d-wave symmetry of the matrix 
element I'k — 4i^7^ is imposed by the hopping geometry: 
a transfer integral v from x^-y^ type Zhang-Rice orbital 
to the neighboring pz states of axial symmetry must have 
different signs along x and y directions. 

Physically, we associate ea with the energy-level sep- 
aration of holes residing on in-plane Op and apical 
Oa sites, as shown schematically in Fig. 1(b). A mag- 
nitude of eA can then be estimated from the Madelung 
potentials on Op and Oa using the structural data 
of Ref. d. To calculate variations of ca caused by super- 
modulation, we use the displacement pattern. Fig. 1(a), 
inferred from the structural data. This way, we quantify 
EA in terms of the distance d from the top Cu02 layer 
(relevant for STM) to Oa above it, making thereby a link 
between the model and supermodulation. 

Below, we use a representative hopping parameters 
t=0.4 eV, t'/t=0.3, t" = t' 12, and vjt = 0.35. The d- 
wave momentum dependence of renormalizes the t', t" 
values but we compensate this numerically by adding a 
counter-term oc v'^ /ex and keep the actual values of t' , t" 
invariant against the Oa shifts (as found in the band 
structure calculations [4]). The covalency effect is then 
entirely due to the spectral weight shifts between x'^-y'^ 
and axial orbitals. 

Calculating the expectation value (c_fc|Cfc|) in the 
model (1), we find the following gap equation: 



1 



— — tanh — — 
2E± 2T 



(2) 



The quasiparticle energies are given by E± 



2iy2±i?2]i/2^ ^here 
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FIG. 1: (a) Local atomic displacements along the c axis 
due to the supermodulation (exaggerated). The Bi-OA bond 
length does not change much but other atoms (Cu, Op, Sr, 
and Ca) are shifted such that the Cu-Oa bond is most affected 
[^. The direction (up/down) and amplitude of the shift 5d 
modulates along the a axis with a period ~26 A, resulting 
in a (5ci ~ ±6% variation of the Cu-Oa bond length 8]. (b) 
Energy level scheme of Cu, Op, and Oa. (c) Schematic pic- 
ture of the screening effect: Once a charge is moved from one 
Cu-ion to another (at the energy cost U, initially) , the apical 
oxygen orbitals experience an electric field (arrows) of excited 
charges and are polarized. An energy gain from the polariza- 
tion process reduces a virtual charge excitation energy U. 



parts where Z± 
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Eq. is composed of two 
A)/iE± — E'i) represent the 



spectral weights of the pairing-active x^-y^ orbital on the 



two bands E±. We note that Z_ (which is the most 
relevant one) is reduced from its bare value {—1 at v = 0) 
due to the orbital mixing. This has the effect of reducing 
effective value of g in Eq. 

First, we regard as a free parameter and consider 
how its variation affects the gap. Solving Eq. ^ at 
g/At — 0.9 and T = 0, we obtained a sizable variation of 
the gap as a function of ea as shown in Fig. 2(a). This 
is due to the covalency effect that reduces the spectral 
weight of the x'^-y^ states near the Fermi-level, by trans- 
ferring it to the higher energy apical states. An amount 
of hole transferred into the apical level is rather small, 
of the order of several percent, see Fig. 2(b). However, 
it may be observed in the nuclear quadrupole resonance 
(NQR) which is sensitive to the hole concentration [22]. 

Next, we consider how ea is shifted by the struc- 
tural modulation. From the Madelung potential calcu- 
lation, we obtained a linear relation £A(d)/i — (ea/OII ~ 
a {d/do - 1)], with £aA=3.85 B and a=1.6. do - 2.4 
A is an average Cu-Oa distance [8|]. The inset of Fig. 2(a) 
shows that the energy level separation ea increases as the 
apical site Oa comes closer to the Cu-ion. Consequently, 
the pairing gap is also increased [see Fig. 2(a)], since the 
strength of the hybridization is reduced. 
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FIG. 2: (a) The pairing gap A as a function of the energy- 
level separation ea between the apical orbital a and the x^- 
band. Arrow indicates ea at d = do. Inset: ea vs. 
Cu-Oa distance d. The data are fitted by a linear relation 
eA{d)/t = 3.85[1 - 1.6(d/do - 1)1- (b) Distribution of doped 
holes among the x^-y^ and the a bands, denoted by pc and pa, 
respectively. Total density of holes is fixed to pc +pa=0.15. 



Having obtained the relations A vs eA and eA vs d, we 
are now in position to show the gap A variations as a 
function of d directly. The result is presented in Fig. 3 
by the broken line, which shows that the ±6% change of 
d leads to a sizable variation of A. 

Screening - This effect is based on an observation 
[l3| that anion polarization renormalizes the energy of 
virtual charge excitations. This physics is relevant here 
since the OA-contribution to the screening of in-plane 
interactions should strongly depend on d. 

Let us consider how U and magnetic correlations 
within the Cu02 planes are modified by apical Oa- In the 
[/-excited intermediate state, both an unoccupied and 
doubly occupied Cu sites strongly polarize apical sites 
just above each Cu site, see Fig. 1(c). Apical oxygen ob- 
tains a dipole moment p — aF where a is the polarizabil- 
ity of O^^ ion, and F ~ e/d^ is an electric field on oxygen 
induced by an extra charge (hole or electron) on Cu-site 
which is created by [/-excitation. (We ignored dipole- 
field corrections to F from further located ions for sim- 
plicity) . Energy gain due to this polarization process in a 
virtual state reduces the [/-excitation energy to Ucsid) = 
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FIG. 3: The gap A as a function of d, where Ao=A(do). 
The broken line shows the covalency effect only, i.e., the gap 
variation is solely due to modulation of eA(d) at fixed g/it = 
0.9. The solid line includes also the screening effect through 
the d-dependence of the coupling constant g{d) (see text). In 
Bi-2212, Cu-Oa distance varies within the range indicated. 



U~Epoi{d), where Epoi{d) = 2(pF/2) = aF'^ ~ a{e/d'^f 
is an interaction energy between the induced dipole mo- 
ments p on Oa and an excited hole (electron) at the 
unoccupied (doubly-occupied) Cu-sites. Since 0~^ ion 
has large polarizability, a ~2 A'^ [2^, energy Epo\{d) is 
sizable: for the average Cu-Oa distance do — 2.4 A we es- 
timate Eq = i?poi('io) —0.9 eV. It is important to realize 
that, due to the strong d-dependence of Epo\{d) oc 
the effective repulsion UeS becomes highly sensitive to 
the Cu-Oa distance: 5Uos{d) = —SEpoi{d) — AEo Sd/do. 
An immediate consequence of this observation is that the 
antiferromagnetic exchange interaction J ~ At'^/UeS be- 
tween the Cu-spins obtains the same strong modulation 
as a function of d: 6J{d)/Jo = — (4i?o/C/o)(d/do - 1), 
where Uq = [/cff(do) and Jq = J{do). With the above 
estimate for Eq and using a representative value Uq ~7 
eV, we find that SJ{d)/Jo = -/3(d/do - 1) with /? ~ 1/2. 
Thus, in a modulated structure of Bi-2212, the J value 
strongly increases as the Cu-Oa distance d decreases and 
vice versa, i.e., it shows the same anticorrelation effect 
with d as the pairing gap A does in the experiment. 

Now, we assume that the superexchange driven mag- 
netic correlations are essential for the pairing in cuprates 
as widely believed. Indeed, e.g., in the t — J model, J 
plays a role of pairing potential same as g in Eq. It 
is then natural to consider that the pairing potential g in 
Eq. ^ is modulated in the same functional form as J . 
We therefore implement a relation g{d) = ^[1 — (3{d/do — 
1)] with /3=l/4 [i^l and calculate the gap values from 
Eq. (21) for different d. The obtained gap modulation is 
presented in Fig. 3 by the solid line. A combined action 
of the covalency and the screening effects can be summa- 
rized by a relation 5A/Aq w —A ■ Sd/do, with A ~ 1.6. 

For a comparison of this result with experiment [lo| . 
we notice that the measured gap A(r) is in fact de- 
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termined by a "coarse-grained" value d{r) « (d(r))j of 
the actual Cu-Oa distances. Such a coarse-graining of 
Sd oc cos(27rrQ/A) gives dd{r) « / • 5d{r), i.e., the 
modulation "seen" by Cooper pairs is reduced by a fac- 
tor / that depends on the ratio of the coherence length f 
and the supermodulation period A ~26 A [25|. At ^ «20 
A, we find / ~ 0.43. The above relation i5A vs Sd reads 
then as SA/Aq k —A ■ Sd/do, where A « 0.7. This gives 
~9% total variation in A due to ±6.25% modulation of 
d, just as observed by Slezak et al. 

In a broader context, it should be emphasized that 
while we are concerned here with the variations of api- 
cal oxygen position within a given structure of a given 
material, the effects discussed - spectral weight transfer, 
and screening of effective U - are generic and relevant for 
the gap and Tc variations among different cuprate fami- 
lies. However, many other things must be kept in mind 
when we compare different cuprates. In particular, apical 
oxygens may have a negative impact on Tc by communi- 
cating a destructive effect of out-of-plane disorder to the 
Cu02 planes (1, 01 . Anticorrelation between the energy 
gap and the distance to the apical oxygen observed by 
Slezak et al. in Bi-2212 implies that the effects we dis- 
cussed here overcome the disorder related physics (which 
would result in a trend opposite to what observed). A 
key question is then how this competition is resolved in 
different cuprate families. To address this issue and bet- 
ter understand the trends in cuprates, our model has 
to be implemented by out-of-plane disorder effects. 

Finally, we argued that effective U hence J values are 
rcnormalized by Oa and thus they become sensitive to 
the OA-position. This inevitably turns the exchange in- 
teraction in Bi-2212 into inhomogeneous one in space. In 
other words, we expect that the strength of local spin 
correlations follow the lattice supermodulation. The re- 
sulting broad distribution of relaxation times could pos- 
sibly be tested by the NMR/NQR experiments. In fact, 
the recent neutron scattering work has revealed an intrin- 
sic broadening of the spin excitations in Bi-2212 [2^, an 
observation that seems natural in a light of our picture. 

To conclude, we discussed the physical origin of the re- 
lationship between the pairing energy gap and the atomic 
displacements in the supermodulated structure of Bi- 
2212. A covalent mixing of the x^-y^ orbital with apical 
p-level, and a screening of effective U values via the polar- 
ization of apical oxygens are found to act cooperatively 
and modulate the pairing correlation as a function of the 
Cu-Oa distance. This leads to the spatial variations of 
the energy gap as observed in the experiment. 
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